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Abstract 

We construct a scaling function /(v^') for inclusive electron scattering from nuclei within the 
Coherent Density Fluctuation Model (CDFM). The latter is a natural extension to finite nuclei 
of the Relativistic Fermi Gas (RFG) model within which the scaling variable ip' was introduced 
by Donnelly and collaborators. The calculations show that the high-momentum components of 
the nucleon momentum distribution in the CDFM and their similarity for different nuclei lead to 
quantitative description of the superscaling in nuclei. The results are in good agreement with the 
experimental data for different transfer momenta showing superscaling for negative values of , 
including those smaller than -1. 

PACS numbers: 25.30.Fj, 21.60.-n, 21.10.Ft, 24.10.Jv, 21.65.+f 
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I. INTRODUCTION 



The ?/-scaling in the inclusive scattering of high-energy electrons from nuclei has been 
actively studied in the last two decades (e.g. 

mm 

and references therein) following 
the idea from It has been shown both theoretically and experimentally that a properly 
defined function (scaling function) depends only on a single variable y, the latter itself 
being a function of the transferred momentum q and energy uj {y = y{q,uj)). It has been 
realized that at momenta |q| > 500 MeV/c and energies u at or below the quasielastic peak 
position a nucleon is ejected in a "quasifree" way almost without the effects of the strong 
interaction. This scaling is usually called scaling of the first kind. It was shown that the 
scaling function is sensitive to the high-momentum components of the spectral function and 
nucleon momentum distribution. Thus its knowledge can provide important information 
about the dynamical ground-state properties, as well as about the reaction mechanism. 

The comparison of the scaling functions of various nuclei with mass number A> 4 led to 
the conclusion that these functions are the same [3, O]- This behaviour is called scaling of 
the second kind which, together with scaling of the first kind, leads to superscaling. These 
studies followed the theoretical concept of the superscaling introduced in P, |l^ considering 
the properties of the relativistic Fermi gas (RFG) model. The analyses of a large body 
of inclusive scattering data for nuclei from A=A to y4=238 in ^ demonstrated that the 
data in the low-cj side of the quasielastic peak exhibit superscaling behaviour: the scaling 
functions are both independent on the momentum transfer and on the mass number. In 
these analyses the Fermi momentum for the RFG was used as a physical scale to define the 
proper scaling variable tp' for each nucleus. An example of the superscaling behaviour of 
the inclusive electron-scattering data for q ^ 1000 MeV/c and for the ^He, ^^C, ^^Al, ^^Fe 



and ^^''Au nuclei is given in Fig^ 
conclusion has been drawn in h 



(the data are taken from Fig. 5 in [8]). An important 
I that this universality is not restricted to the region 
of the quasielastic peak (lip'l < 1) and that the superscaling extends to larger values of 
and thus, to the high- momentum components of the nucleon momentum distributions 
in nuclei. The existence of high-momentum components, and their similarity for different 
nuclei, is known to be due to the short-range and tensor nucleon-nucleon correlations (see, 



e.g. 



and references therein). 



An extended study of scaling of the first and second kinds for inclusive electron scattering 
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from nuclei with emphasis on the transverse response in the region above the quasielastic 



peak was performed in 



111 ]. Approximate scahng of the second kind was observed and its 



modest breaking was supposed to be due to an inelastic version of the usual scaling variable. 



In [1^ a unified relativistic approach used in the case of quasielastic kinematics was applied 
to the analysis of highly-inelastic electron scattering. The complete inelastic spectrum was 
considered using the inelastic RFG model and its phenomenological extension based on 
direct fits to data, investigating the second-kind scaling behaviour as well. 

As emphasized in jsl , the actual dynamical physics content of the phenomenon considered 
is more complex than that provided by the RFG framework. In particular, as noticed there, 
the extension of the superscaling property to large negative values of ip' < —1) is not 
predicted by the RFG model. This is seen in Fig. Q where we also give a curve of the 
calculated RFG scaling function which is equal to zero for ip' < —1. Thus, it is worth 
considering the superscaling in theoretical approaches which go beyond the RFG model. 
One of them is the Coherent Density Fluctuation Model (CDFM) (e.g. H Q Q, Q) 
which gives a natural extension of the Fermi-gas case to realistic finite nuclear systems. The 
main aim of the present work is to see to what extent superscaling can be explained using 
the CDFM. The theoretical scheme is given in Section II, while the results and the discussion 
are presented in Section III. The conclusions and the final remarks are given in Section IV. 

II. THE THEORETICAL SCHEME 
A. Basic relationships of the CDFM 

The CDFM was suggested and developed 14, 15, 16, [l7| as a model for studying 



characteristics of nuclear structure and nuclear reactions based on the local density and 
momentum distributions as basic variables of the theory and using the essential results 
of the infinite nuclear matter theory. The model is related to the delta-function limit of 
the generator coordinate method (GCM) |18]. In the latter the total many-particle wave 
function \E'({rj}) of a system of A nucleons is written in a form of a linear combination: 

^(ri, ...,rA) = J F{xi,X2, r^; Xi, 0:2, ...)dxidx2..., (1) 

where the generating function $({rj}; Xi,X2, ...) depends on the radius-vectors of the nucle- 
ons {rj} (spin and isospin variables are implied) and on the generator coordinates Xi,X2i .... 
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FIG. 1: Superscaling behaviour of inclusive electron-scattering. The grey area represents experi- 
mental data 0] for "^He, ^^C, 27 ai and ^^^Au ai q = 1000 MeV/c. The solid line is the RFG scaling 
function calculated using Eq. (EH) with /ciT = 1.191 fm"^ from Ref. 

The function $ is usually chosen to be a Slater determinant built up from single-particle 
wave functions corresponding to a given construction potential parametrized by xi^X2i---- 
The weight function F(xi,X2,...) can be determined using the variational principle as a 
solution of the Hill- Wheeler integral equation: 



j x') - EI{x, x')]F{x')dx' = 0. 



(2) 



In Q the overlap kernel I{x, x') and the energy kernel T-C{x, x') have the following forms, 
respectively: 

J(x,a;') = (^'({r,},a;)|$({r,},x')), (3) 
n{x,x') = (<l>({r,;},x)|^|$({r,},x')), (4) 



where i=l,2,...,y4, H is the Hamiltonian of the system and x denotes a set of xi , x^-^ .... For 
many-fermion systems the kernels 1{x, x') and 7i(x, x') peak strongly at x ~ x' 



Xi^^ 

c' [igiy 



can be written in the form: 



X(x, x') ^ X(x, x)Q{x — x'), 
'H{x, x') ~ H(a;, x)Q{x — x'), 



and 

(5) 
(6) 



where Q is peaked at x ~ x'. It was shown in that the following delta-function approx- 
imation for the kernels is valid in the GCM for the case of many-fermion systems 

I{x,x') 6{x-x'), (7) 
n{x, x') -> -^r^S"{x -x') + V (^^] 6{x - x') (8) 



and that it leads to the Schrodinger-type of equation [18|, |l9|, |2l| with an effective mass 
dependent on the generator coordinate (see also 

In the following we use for simplicity only one generator coordinate. If the trial wave 
function \l/({ri}) in is normalized to the mass number A and the weight function is 
determined under the condition 



00 

2, 







\F{x)\^dx = 1, (9) 

then the delta-function approximation ((Tj) leads to the relationship: 

J r ({r^}, x')$({r,}, x)dri...drA = AS{x - x'). (10) 

Taking into account Eqs. (jSj) and (fTO|) it was suggested in the CDFM that a delta-function 
approximation leading to that of Eq. holds in the case of many-fermion systems [14-16]: 

J $*(r,r2, ...,rA,a;')$(r',r2, ...,rA,x)dr2...drA = Px,xi^,r')6{x - x'). (11) 

In Eq. ()11|) Px^x{t^,t^') is the one-body density matrix corresponding to the wave function 
^({rj},^) which can be formally written as 

Px,x(r,r') = px{r,r') 
A r 

= J ^*{r,r2,...,rA,x)^{r',r2,...,rA,x)dr2...drA. (12) 

One can see that the integration of Eq. (fTT|) (at r'=r) over r using Eq. ()12|) leads to Eq. 
(fTUI) which is the delta-function limit for the overlap kernel in the GCM. 

In the CDFM the generating function <l>({rj},x) describes a system corresponding to a 
piece of nuclear matter with a one-body density matrix (ODM) of the form 

/ /X o / ^ j^ikF{x)\r -r'\) ( |r + r^| ^ 
P^^'^ ' ) = '^"^"^ {Mx)\r-r'\) Q - ) 

and uniform density 

Px{r) = po{x)e{x-\r\), (14) 



where 



Po[x) 



3A 



(15) 



and the generator coordinate x is the radius of a sphere containing all A nucleons in it. 
In Eq. (fT^ ji is the first-order spherical Bessel function and 



1/3 



-po[x) 



a 

X 



with a 



( 



1.52^1/3 



(16) 



is the Fermi momentum of such a piece of nuclear matter. Using Eqs. (^J and the 
ODM of the systeni in the CDFM can be obtained as a superposition of the ODM's from 
Eq. (USD 



14 



15 



16 



17|: 



p(r,r')= / dx\F{x)\'p,{r,v'). (17) 



The Wigner distribution function which corresponds to the ODM from Eq. (|T7jl is: 

dx\F{x)\'W,{r,k), (18) 



where 



W{r,k) 



(2vr)' 



re(x- |r|)e(fcF(x) - |k|). 



(19) 



Then the density p(r) and the momentum distributions n(k) in the CDFM are expressed 
by means of the same weight function |F(x)p: 



/roo Q /I 

dkW{r,k) = J^ dx\F{x)\^^—^e{x-\r\) 



(20) 



and 



n(k) 



drW{r,k) 



(2vr) 



dx\F{x) 



Atxx^' 



e{kF{x) - \k\) 
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(27r)3 

both normalized to the mass number: 



/ dx\F{x)\^-Tix'' 
Jo 3 



J p{r)dr = A; j n{k)dk = A. 



(21) 



(22) 



Different paths can be followed to find the function F{x). Here, instead of solving the 
differential equation from the delta-function approximation to the Hill- Wheeler integral 
_we adopt a convenient approach to the weight function F{x) proposed in 
The function F{x) is obtained by means of a known density distribution 
p(r) for a given nucleus [from Eq. ^K^]: 




\F{x) 



1 dp{r) 
Po{x) dr 



(at dp{r)/dr < 0). 



(23) 
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B. The scaling function in the CDFM 



The scaling function in the RFG model expressed by the variable if)' has the form 



2 + r]^p~2Jl + r]j. 



(24) 



where rjp = kp/rriN, mjq being the nucleon mass. 

As shown in , the relationship between ifj' and the usual variable, in the approximation 
for the mass of the residual nucleus Mj^_^ ^ oo, is given by the expression 



y_ 



1 



(25) 



where k = q/2mjq. 

Our basic assumption within the CDFM is that the scaling function for a finite nucleus 
f{ip') can be defined by means of the weight function |F(a;)p, weighting the scaling function 
for the RFG at given x (i.e. for a given density Pq{x) (jl5p and Fermi momentum 
Thus the scaling function f{ip') in the CDFM will be an infinite superposition of the RFG 
scaling functions f{ip'{x)). 

Let us introduce the notation 



^ _ 1 / ^ 1 



(26) 



Then one can write from Eqs. ()25p and (j^Bj) . neglecting the scaling variable ip'^iy) 

corresponding to the relativistic Fermi gas with the density po{x) (jlSj) and Fermi momentum 
kpix) in the form: 



p{y) ^ p{y)x 

kp{x) a 



(27) 



where for the cases of interest 



p{y) 



y{^ + cy), y>0 
-\y\il-c\y\), y<0, < l/(2c). 



(28) 



For further use it is more convenient to introduce the notation 

kp p{y) kp 



kp{x) kp kp{x 



(29) 
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Using the 0-function in Eq. ()24|) . the weighted scahng function for a finite nucleus can be 
presented by the integral 



dx\F{x)\'- 



\ a 



X 



a J 



2 + 



/ a 



2W1 + 



/ a 



yxniN 



(30) 



where the Fermi momentum kp will not be a fitting parameter (as it is in the RFG model) 
for the different nuclei, but will be also calculated consistently in the CDFM, 



kp = dxkF{x)\F{x)\'^ = a / dx—\J-'{x) 
Jo Jo X 



(31) 



for each nucleus, with a given by Eq. (fTB|) . As can be seen from Eqs. (IHUj). and (^Hj) 
in our approach the scaling function f{ip') is symmetric at the change of ip' to -ijj' up to 
l^'l = {4ckp)-\ 

The scaling function f{ip') has been calculated using Eq. ()30|) by means of the weight 
function |F(x)p determined from its relationship to the density distribution p(r) [Eq. (PHj) ]. 
For the latter we used those obtained from experimental data on electron scattering from 
nuclei and muonic atoms. 



III. RESULTS OF CALCULATIONS AND DISCUSSION 

We calculated the scaling function f{ip') (jHUj) for various nuclei and transfer momenta. 
A symmetrized diffused Fermi density distribution has been used for ^He and ^^C [22] and 
a diffused Fermi distribution for the heavier nuclei. The values of the half-radius R and 
diffuseness parameter b are given in Table H] together with the results for the CDFM Fermi 
momentum kp fl31|) . 

The results for the scaling function are compared with the experimental data from 0, 
which are given in our figures by a grey area. In Fig. |21 are presented the results for the 
scaling function in the CDFM for g=1000 MeV/c and for ^He, ^^c, 27A1 and ^^^Au. The 
values of the parameters R and b for ''He and ^^C (given in Table 1) lead to ch arg e rms 
radii 1.71 fm and 2.47 fm, respectively, which coincide with the experimental ones [2^. The 



values of R and b for ^''Al are taken from 



23|. The results of the CDFM scaling function 



(solid lines) are compared with the RFG predictions (dotted lines). In the RGF model, due 



TABLE I: Values of the parameters R and h (in fm) used in the calculations and the results for kp 
(in ixnr^) obtained in the CDFM. 



Nuclei 


R h kp 


^He 


1.710 0.290 1.201 


12 c 


2.470 0.420 1.200 


27A1 [23j 


3.070 0.519 1.267 


^^Fe [23] 


4.111 0.558 1.270 


is^u [24] 


6.419 0.449 1.335 



to the B-function in Eq. ()24|). f{ip') = for V'' < — 1. As can be seen, the CDFM results 
give a good agreement with the data for the interval for ip' from till ip' < —1 for all nuclei 
(including ^^Fe which is not shown). The only exception was observed for ^^^Au using the 
values of the parameters R = 6.419 fm and b = 0.449 fm given in for which the result 
is shown in Fig. |21by dashed line. 

At this point we would like to note that, generally, the weight function |F(x)p which is 
related to the density distribution [Eqs. (plj) and (f^ ]. is also related to the momentum 
distribution [Eq. (j2T| ]. This connects through Eq. ()30|1 the scaling function f{ip') with 
n{k). The deviation of fiip') from the data in the case of ^^'^Au {b = 0.449) is due to the 
particular A-dependence of n{k) in the present approach [Eqs. (pT|l . pHjl ]. 

To understand the origin of the agreement for light and medium nuclei, and of the dis- 
crepancy for ^^^Au with b = 0.449 fm, we show in Fig. El the momentum distributions n{k) 
(EH) (Fig- Efa)) and the weight functions |F(x)p (Fig. Ofb)). As seen in Fig. Efa), for the 
^^C and "^^Ca the CDFM momentum distributions depend weakly on A and have similar 
high-momentum {k > 1 fm~^) tails. Inspection of Eqs. (f7T|l and (jHn|) (which involve the 
same weighting function |F(x)p) shows that the interval of interest for ip' (—2 < ip' < 0) 
corresponds to the k interval < A; < —2.5 fm~^ in n{k). Thus the similar high-momentum 
tails of n{k) in the light (^^C) and medium (^°Ca) nuclei lead to the similar CDFM scaling 
curves which are in agreement with the superscaling data. 

The deviation between the CDFM scaling function and data in Fig. |2lfor ^^^Au {b = 0.449 
fm) takes place at ip' < —1 and corresponds to values of the momentum A; > 3 fm~^. The 
CDFM high-momentum tail of ^^''Au for b = 0.449 fm is much smaller than that of light 
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FIG. 2: Results for the scaling function in the CDFM (solid line) calculated using Eqs. and 
(ISU) atq= 1000 MeV/c and for ^He, ^^C, ^'^Al and ^^^Au (with b = 1.0 fm for the latter) compared 
with the data (grey area) from j^. The dotted line is the RFG result using Eq. (|24() . The dashed 
line in the case of ^^''Au corresponds to the CDFM result with b = 0.449 fm. 

and medium nuclei (see Fig. Ofa)) causing the above mentioned deviation. In Fig. ^h) we 
give also the weight function |F(x)p calculated in the CDFM for ^He, ^^C, ^°Ca and ^^^Au 
using Eq. (|23|) and the corresponding density distributions mentioned above. As can be 
seen |F(a;)p is a one-peak function which for ^He, ^^C and ^''Ca follows a particular trend 
with increasing A. The strength of the peak decreases (consequently its width increases) 
with increasing A, and the peak is displaced to higher x values. This trend is broken by the 
behaviour of |F(a;)p for ^^^Au when we take b = 0.449 fm. 

To improve the A-dependence of the momentum distribution in the CDFM for the heav- 
iest nucleus, we use a procedure that may be somewhat artificial but which is useful to show 
the role of the obtained, more realistic, new high-momentum components of n{k) on the 
scaling function. This can be achieved by taking an effective larger value of the parameter b 
in the diffuse Fermi density distribution which is used to obtain the weight function |F(a;)p 
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FIG. 3: (a) Nucleoli momentum distribution n(k) calculated in the CDFM using Eqs. (|2H) and 
((21 for 12c, ^°Ca and ^^'^Au (for the latter with b = 0.449 fm and b = 1.0 fm); (b) The weight 
function \F{x)\^ of the CDFM calculated using Eq. (El for "^He, ^^c, ^OCa and i^'^Au (for the 
latter with b = 0.449 fm and 6 = 1.0 fm). 

and, hence, the scahng function f{ip'). We take the value b = 1.0 fm, for which the high- 
momentum components of n{k) in ^^''Au are similar to those in light and medium nuclei. 
This can be seen in Fig. Efa) for ^^^Au (solid line). The change of the behaviour of the 
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weight function which leads to this shape of n{k) for ^^^Au can be seen in Fig. Efb) 

(sohd hne). In this case the function follows the trend previously observed for light 

and medium nuclei of decreasing strength (increasing width) with increasing A. The values 
of for < a; < 5.2 fm determine the behaviour of f{ip') for —1.5 < ip' < —1. As 

can be seen from Fig. El^b) in this region is quite different for b = 0.449 fm and for 

b = 1.0 fm. The results of the calculations of the scaling function for ^^^Au with b = 1.0 fm 
are presented in Fig. |21by a solid line and, as can be seen, they are in good agreement with 
the data. This confirms the view that the behaviour of the scaling function is related to the 
properties of the momentum distribution at large values of k {k > 1.5 fm~^), and that the 
similarity of the high-momentum tails of n{k) leads to the scaling of second kind. 

Here we would like to note that the use of an effective value of b for ^^^Au can be merely 
seen as an artificial procedure to improve the A-dependence of n{k) for the heaviest nucleus. 
This shows what would be the results of the CDFM for the scaling function when the high- 
momentum tails of n{k) are realistic, even for the heaviest nucleus, and are similar to those 
of light and medium nuclei. We do not imply that the actual diffuseness of the density 
distribution of ^^^Au should be that large. However, it is also worth pointing out that all 
the nucleons may contribute to the scaling function for the transverse electron scattering 
and that the diffuseness of the mass density for a nucleus like ^^'^Au may be larger than that 
of the charge density. 

The results for the scaling function in the CDFM in the case of g=1650 MeV/c are given 
in Fig. 13 for ^He and ^^''Au (the latter with improved high-momentum tail of n{k)) and 
compared with the experimental data taken from Fig. 6 of The curves for ^^C and ^^Fe 
are not given since they are similar, in agreement with the data. 

The result for the scaling function in the CDFM in the case of g=500 MeV/c for ^^C is 
given in Fig. |S| This result is also in good agreement with the experimental data for q in 
the interval from 500 to 600 MeV/c given in Fig. 8 of js!]. 

One can see in Figs. El and E] that the CDFM results tend to overestimate the data in the 
interval — 1 < V^' < —0.5. We note that the origin of this is related to the predictions of the 
RFC model in the same region, as can be seen in Fig. 
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FIG. 4: Results of the CDFM for the superscahng functions of ^He (soHd hne) and ^^^Au (dashed 
hne) aX q = 1650 MeV/c compared with the experimental data (grey area) from Q]. 

IV. CONCLUSIONS AND FINAL REMARKS 

The results of the present work can be summarized as follows: 

i) We propose an extension of the RFG model to calculate the scaling function f{ip') in 
finite nuclei within the Coherent Density Fluctuation Model (CDFM). In this model f{ip') 
is a weighted superposition of scaling functions for relativistic Fermi gases with different 
densities. The weight function is calculated using the known charge density distributions in 
nuclei. 

ii) We calculate the scaling function f{il)') for inclusive electron scattering for ^He, ^^C, 
^^Al, ^^Fe and ^^^Au nuclei and for various values of the transfer momentum |q|=1650, 1000 
and 500 MeV/c. The results agree with the available experimental data at different trans- 
ferred momenta, and energies below the quasielastic peak position, showing superscahng for 
negative values of ip' including also those smaller than -1. This is an improvement over the 
RFG model predictions where the scaling function becomes abruptly zero beyond = —1. 

iii) The sensitivity of the scaling function to the high-momentum components of the 
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FIG. 5: Results of the CDFM for the superscahng function of ^^C at q = 500 MeV/c (soHd hne) 
compared with the experimental data (grey area) for q in the interval from 500 to 600 MeV/c from 

momentum distribution is analyzed in detail, especially so on the example of the ^^''Au 
nucleus. 

iv) The scaling function obtained is symmetrical around = up to = {Ackp)'^- 
It would also be interesting to search for models predicting an assymetrical superscaling 
function /{"ip') as the phenomeno logical one obtained by ^25] and discussed in Ref. [l 2| . 

It is shown in our work that the superscaling in nuclei can be explained quantitatively 
on the basis of the similar behaviour of the high-momentum components of n{k) in light, 
medium and heavy nuclei which is known to be due to the short-range and tensor correlations 
in nuclei. This suggests an alternative path for defining the weight function F{x) within the 
generator coordinate method: a path in which F{x) is built up from a phenomenological or 
a theoretical momentum distribution. 
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